Wigner Crystals of 229 Th for Optical Excitation of the Nuclear Isomer 
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We have produced laser-cooled Wigner crystals of 229 Th 3+ in a linear Paul trap. The magnetic 
dipole (A) and electric quadrupole (B) hyperfine constants for four low-lying electronic levels and the 
relative isotope shifts with respect to 232 Th 3+ for three low- lying optical transitions are measured. 
Using the hyperfine B constants in conjunction with prior atomic structure calculations, a new value 
of the spectroscopic nuclear electric quadrupole moment Q = 3.11(16) eb is deduced. These results 
are a step towards optical excitation of the low- lying isomer level in the 229 Th nucleus. 



Atomic nuclei usually possess excitation energies in the 
keV to MeV range. A notable exception is the 229 Th nu- 
cleus where the energy splitting of the ground state dou- 
blet is less than 10 eV [ESI- This fortuitous coincidence 
offers prospects for a better clock and for a test of tem- 
poral variation of fundamental constants at an unprece- 
dented level of precision [5] [6] . Here we report laser cool- 
ing and crystallization of 229 Th 3+ in a linear Paul trap. 
The monovalent character of the Th 3+ ion is favorable 
for the isomer level search based on an electron-bridge 
process [7]. Once the isomer energy value is measured, 
either in our system or using one of other approaches 
[8] [9], a single trapped, cold 229 Th 3+ ion is expected to 
be an ideal system to take advantage of this remarkable 
nucleus. 

Because to the intrinsic narrowness of the nuclear tran- 
sition and the superb isolation of the nucleus from exter- 
nal fields, 229 Th may form the basis of a next-generation 
optical clock [5 . To this end, a convenient set of elec- 
tronic transitions in the triply charged form of the Th 
atom may be utilized for optical preparation and readout 
of the nuclear ground and isomer manifolds. The clock 
levels may be associated with the metastable 7S elec- 
tronic levels of the ground and isomer manifolds, while 
allowed electric dipole transitions from the F states are 
used for laser cooling of the ion and for optical pumping 
preparation of the clock, Fig. [I] 

Importantly, the ground-isomer nuclear transition may 
be used to search for temporal variation of fundamental 
constants by monitoring the ratio of the frequencies of 
the nuclear transition and that of an atomic (electronic) 
clock transition, e.g. in Al + or Sr. Flambaum et al 
have predicted that such a frequency ratio is likely to be 
supremely sensitive to relative variation of the strong in- 
teraction parameter and the fine structure constant, with 
enhancement factor K possibly reaching as high as 10 6 
[6] [lOl [11] . Although the complexity of this nuclear sys- 
tem makes a theoretical evaluation of K difficult p~2j [13] , 
optical and microwave spectroscopy of the ground and 
isomer level manifolds should provide a reliable empiri- 
cal determination of the enhancement factor with little 
or no dependence on the nuclear model assumed [13] . 

In this Letter, we report an important step toward 
these goals: laser cooling, to crystallization, of 229 Th 3+ 



in a linear Paul trap, with the ultracold ions stored for 
a large fraction of an hour. Triply charged 229 Th and 
232 Th ions are created via ablation [a 100 /i J, 5 ns pulse at 
the third harmonic of a yttrium aluminum garnet (YAG) 
laser] of a thorium nitrate source and injected axially 
into a four-rod linear Paul trap, Fig. [2] The design of 
the trap is motivated by the plasma nature of the Th 3+ 
ablation plume and is tailored for enhanced ion loading 
efficiency [T5] , The trap used here is of length 188 mm 
and radius 3.3 mm, with a 10 mm distance between the 
thorium nitrate source and the trap. 

The trap is driven in the balanced mode at 4.3 MHz 
and single-ion radial trap frequencies for both isotopes 
are varied from 10 kHz - 1 MHz and the axial frequen- 
cies between 15 - 55 kHz during experimental sequences. 
To create sufficient axial confinement with convenient dc 
voltages, two tube electrodes, of length 68 mm and 83 
mm, sit coaxially with the rf potential and are separated 




FIG. 1. A diagram of electronic energy levels and electric- 
dipole transitions including both nuclear ground and excited 
isomeric manifolds. Direct nuclear magnetic dipole (Ml) de- 
cay channels and electron-bridge pathways suitable for the 
isomer search are shown. Optical transition wavelengths are 
in nm and the integers near atomic levels indicate principle 
quantum numbers. 
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FIG. 2. Essential elements of the experimental setup. See 
text for discussion. 




FIG. 3. Various optical transitions between hyperfine levels in 
229 Th 3+ employed during laser cooling. See text for details. 



by 22 mm at trap center. The electrodes are biased at 
3-50 V. A tapered geometry of the tube inner diameters 
creates a potential gradient along the trap axis, inducing 
ion accumulation between the tubes. For ion loading, a 
6 mm long ring electrode biased at 250 V sits at each 
extreme end of the trap. The ring electrode voltages are 
pulsed down to ground and back in sync with ion abla- 
tion, allowing Th 3+ entrance to the trap while blocking 
much of the slower moving Th 2+ and Th + . The mass- 
dependent stability of rf trapping is also utilized to avoid 
trapping pollutant ions and to remove pollutants from 
the trapped sample [T5] . 

Once trapped, laser cooling is sufficient to bring the 
ions from tens of thousands of Kelvin to crystallization, 
though overall yield is consistently < 100 ions. To in- 
crease this yield by 1 - 2 orders of magnitude, helium 
buffer gas of impurity < 100 ppb is introduced into the 
vacuum system at a pressure of 10 ~ 5 Torr, providing ini- 
tial cooling to about room temperature. The buffer gas 
is then removed, allowing for laser cooling to take place. 
The vacuum environment returns to a base pressure of 
< 10 -11 Torr within a few minutes thereafter. 

Laser cooling of trapped Th 3+ isotopes proceeds in two 
stages. First, transitions power broadened to several hun- 
dred times their natural linewidths are used for cooling 
ions from tens of thousands to tens of Kelvin. The rather 
long excited state lifetimes of the 6D 3 / 2 an d 6D5/2 levels 
- of order 1 /is [16J - are compensated for by the high laser 
intensities to allow cooling to overcome the heating asso- 
ciated with micromotion. Subsequently, laser intensities 
which power broaden transitions to only several natural 
linewidths are utilized in cooling the ions into the mK 
regime. 

To scatter light from and laser cool 229 Th 3+ in the 
steady state, multiple optical fields are required to re- 
pump from the various ground (5F 5 / 2 ) and metastable 
(5F7/2) hyperfine levels. Fig. [3] illustrates relevant en- 
ergy levels of 229 Th 3+ and the scheme used in this work; 
the 1088 nm transitions are used in the high-temperature 
regime and the closed \5F 5 /2,F — 1)> |6£>5/2,^ = 0), 



|5F 7 / 2 ,^ = 1) A system is used in the low-temperature 
regime. In this case, the requirement of multiple finely 
tuned optical frequencies is relaxed, as only repumping 
by 1088 nm and 984 nm fields is needed to correct for 
off-resonant excitation into the 6D5/2 manifold. The 
232 Th 3+ system has no hyperfine structure, and there- 
fore, no such complications arise [7J. 

The 1088 nm transitions of both isotopes are uti- 
lized for cooling the sample down to tens of Kelvin 
using two diode lasers. The 229 Th 3+ light is initially 
passed through a fiber-based electro-optic phase modu- 
lator (EOM), acquiring rf sidebands in the 0.1-1 GHz 
range. In total, four 2 mm diameter beams, oriented in 
a common plane, converge on the ions at 45° angles to 
the trap axis. The 232 Th laser and the 229 Th laser carrier 
and rf sidebands are tuned such that spectral components 
are located 20-80 MHz below relevant transition frequen- 
cies. Each spectral component addressing 229 Th 3+ car- 
ries about 1 mW of power per beam and the 232 Th 3+ 
field carries about 20 mW per beam. The lower 229 Th 
power is a result of EOM optical insertion loss and the 
distribution of laser power over multiple rf sidebands. 

For low-temperature cooling, transitions from the 
5F 7 /2 manifolds to the 6D5/2 manifolds are excited with 
light from a frequency-stabilized 984 nm diode laser 
which passes through an EOM, acquiring rf sidebands in 
the frequency range 3-7 GHz. The spectral components 
are placed 3 MHz below the resonant frequencies of the 
232 Th 3+ and relevant 229 Th 3+ transitions. The light is 
then directed to the ions with two 1 mm diameter beams 
counterpropagating along the trap axis. Each beam con- 
tains about 20 /iW of power per spectral component. 
This beam orientation minimizes observed micromotion- 
induced Doppler broadening. A nearly identical system 
at 690 nm is used to excite from the 5F 5 / 2 manifolds to 
the 6D5/2 manifolds. Each 690 nm beam carries about 60 
/iW of power in each relevant spectral component. When 
ion temperature and Doppler width are reduced via ini- 
tial 1088 nm laser cooling, the A system fields become 
relevant, cooling both isotopes to crystallization. 
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During the laser cooling process, 984 nm fluorescence 
is used for ion detection. Light is collected with a high 
numeric aperture (NA = 0.55) aspherical lens and di- 
rected onto a low-noise CCD camera. The collection effi- 
ciency here is about 7%. Fig. [4] shows images of various 
Wigner crystals of 229 Th 3+ and 232 Th 3+ . Images were 
taken with different 690 nm fields present, controlling 
the level of 984 nm fluorescence from each of the two iso- 
topes. The 10.5 GHz isotope shift is much larger than 
the broadened cooling transition widths of a few MHz, 
allowing for unambiguous isotopic identification. The rel- 
ative abundance of the two isotopes (76% 229 Th and 24% 
232 Th) in the thorium nitrate source is reflected in the 
observed composition of the dual-isotope crystals. When 
low-temperature laser cooling is not present for one of the 
two isotopes, an ultracold phase is maintained due to the 
apparent sympathetic cooling. For single ion excitation 
and detection, the most robust geometric configuration 
is the linear ion chain, such as those shown in Fig. |4]3,d. 
Ions in this configuration are observed to survive in the 
trap for a large fraction of an hour. The lifetime is likely 
limited by charge exchange collisions and chemical reac- 
tions with background molecules, processes enhanced by 
the large electron affinity of Th 3+ [17] . 

To determine and manipulate the internal states of ul- 
tracold 229 Th 3+ , optical transition frequencies between 
the multitude of hyperfine levels must be known at the 
level of a few MHz. These frequencies are measured using 
optical pumping and fluorescence techniques via suitably 
engineered spectra at all three wavelengths. To clearly 
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FIG. 4. Trapped 229 Th 3+ and 232 Th 3+ ions, laser cooled to 
crystallization, (a) Three successive images were taken of 
about 200 crystallized Th 3+ ions. The top image was taken 
with both 229 Th and 232 Th 690 nm fields present. The 232 Th 
690 nm field was then removed for the middle image. The 
232 Th 690 nm field was then reintroduced and the 229 Th 690 
nm field simultaneously removed for the bottom image. The 
larger mass- to- charge ratio of 232 Th 3+ causes its radial ac- 
cumulation in the outer shells of the crystal. The short fo- 
cal depth of the imaging system is apparent as most of the 
232 Th 3+ shell is out of focus. The integration time is 1 s for 
all three images, (b), (c) The same imaging protocol used in 
(a) is used on smaller samples. The integration time is 2 s for 
all images, (d) A linear chain of four 229 Th 3+ ions is shown. 
The integration time is 3 s. The scale bar in (a) is 500 fim. 
and in (b)-(d) are 100 /xm. 



resolve all optical transitions, the ions are sympatheti- 
cally cooled by 232 Th 3+ to tens of mK, narrowing the 
transition Doppler widths to a few MHz. A single hyper- 
fine level in the 5F 5 / 2 or 5F 7 / 2 manifold is left dark so 
that the ions are optically pumped into it. The 984 nm 
fluorescence is monitored on a camera while a probe field 
frequency is scanned. When this field is resonant with the 
targeted transition, the dark level is depopulated and a 
concomitant fluorescent signal is observed. At maximal 
fluorescence, the value of the probe frequency is recorded 
with respect to the corresponding 232 Th 3+ transition. 

The magnetic dipole (A) and electric quadrupole (B) 
hyperfine coefficients of all four levels and the relative 
isotope shifts of all three transitions are extracted di- 
rectly from the recorded data via a global least-squares 
fit containing expressions of the form [18] 

_ _Ke, , |^ e (^ e + l)-2/(/+l)J e (J e + l) D 

°^ 9 - 2 Ae + 4/(2/ - 1) J e (2J e - 1) 

_^, 4 %K g (K g + l)-2I(I + l)J g (J g + l) 

2 9 41(21 -l)J g (2J g -l) 9 ^ 

where the subscript g(e) indicates the ground (excited) 
state, Ki = Fi(Fi + 1) - /(/ + 1) - Ji(Ji + 1), and A 
represents the relative isotope shift common to all optical 
transitions between a given pair of hyperfine manifolds. 
Results of the fit are presented in Table [TJ 

Overall, about 20% of the uncertainty is statistical for 
each extracted value, while the remaining 80% is due 
to uncertainties in Zeeman splittings combined with op- 
tical pumping effects, laser cooling dynamics, and light 
shifts. By taking the unweighted average ratio of the four 
measured B's to the corresponding calculated electronic 
hyperfine matrix elements of Ref. [14], a value of the 
spectroscopic nuclear electric quadrupole moment Q = 
3.11(16) eb is obtained. Here, the uncertainty is deter- 
mined by the few percent accuracy of the atomic struc- 
ture calculations. No accurate calculations are currently 
available to extract the value of the nuclear magnetic 
dipole moment, although they are possible [19]. Our 
value for Q may be compared with the previous value 
of Q = 4.3 eb ± 20% measured in Th+ [2D]. 

The most recent and precise published measurement 
of the 229 Th isomer energy is 7.6(5) eV [4]. In order to 
span ± 3 a in the search for this nuclear level, direct op- 
tical excitation of the isomer in trapped cold ions may 
not be a viable method, given available UV sources and 
the large energy uncertainty. Instead, the electron-bridge 
(EB) process may be utilized [7]. In this case, hyperfine- 
induced mixing of the ground and the isomer nuclear 
manifolds opens electric-dipole transitions between the 
two. The mixing is expected to be the strongest for the 
5-electronic states, as the electron probability density at 
the nucleus is highest. For example, considering only the 
7Si/2 and 851/2 electronic orbitals in first-order pertur- 
bation theory, 

\7S 1/2 ,m) « |75i /2 ,ra) + \SS 1/2 ,g) 
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TABLE I. Measured 229 Th 3+ hyperfine constants A and B 
and relative isotope shifts from 232 xh 3+ . All units are MHz 
and all uncertainties are one sigma. 



Valence 








Isotope 


Orbital 


A 


B 


Transition 


Shift 


5F 5 /2 


82.2(6) 


2269(6) 


5F 5/2 6D 3/2 


-9856(10) 


5F 7/2 


31.4(7) 


2550(12) 


5F 5/2 6D 5/2 


-10509(7) 


6£> 3 /2 


155.3(12) 


2265(9) 


5F 7/2 6L> 5/2 


-9890(9) 


6A3/2 


-12.6(7) 


2694(7) 







where i^ n £ is the electron-nucleus interaction Hamilto- 
nian and g(m) indicates the nuclear ground (isomer) 
level. The |8Si/2,g) admixture, with expected ampli- 
tude of order 10 -5 [21], couples to the \lP\/2,9) level via 
electric-dipole radiation of frequency (E 7 s^m—^7P 1/2 ,g)/^ 
(see Fig. [I]). This shifts the spectral interrogation re- 
gion from the challenging 130-200 nm range to the man- 
ageable 250 - 800 nm range where high-power coherent 
light sources are available. It should be noted that this 
electron-bridge transition is predicted to be substantially 
stronger than the nuclear Ml transition [2T] . 

In general, the optical search time required to find 
the isomer transition energy can be estimated as T ^ 
Auj/Q 2 , where Q is the resonant optical excitation Rabi 
frequency and HAuj is the energy range over which 
the search is conducted. Considering the |7Pi/2,#) 
1 75i/2, m) electron-bridge, an effective electric-dipole 
moment, (Ieb, is predicted to be ^ 2 x 10 -5 eao 
[21]. For an ion in the \lP\/2,9) level, placed at 
the center of a Gaussian light beam of power P and 
waist r , ft 2 = (4P (Ieb 2 )/^ e cr 2 ^ 2 ) so that T ~ 
(7T6ocr 2 h 2 Auj)/(4Pd EB 2 ). 

The spectral range in which this transition is likely 
to exist (250 nm-800 nm) is completely accessible with 
a mode-locked Ti:Sapp oscillator at the fundamental (1 
W), second (150 mW), and third (20 mW) harmonics, 
and with an optical parametric oscillator in the visi- 
ble spectrum (200 mW). By focusing ultra- fast excita- 
tion pulses from these or similar light sources to a waist 
ro ~ 5 /ira and onto a sing le crystallized 229 Th 3+ ion, we 
estimate that the isomer level can be excited from the 
\7Pi/2,g) level in tens of hours of ion illumination. Alter- 
natively, \7Si/2,g) \7Si/2,m) two-photon excitation 
through the |7Pi/ 2 ,g) level may be a suitable approach 
for the search. It should be noted that the ionization po- 
tential of Th 3+ is 29 eV, making multi-photon ionization 
with these protocols unlikely. 

Once the isomer level is found, a single ion within a lin- 
ear crystallized chain (Fig. [4^,01) can been excited to the 
isomer manifold and its hyperfine structure and isomeric 
level shifts may be accurately measured. This knowledge 
would allow both for unambiguous identification of the 
isomer level and for an empirical determination of the 
isomer transition sensitivity to a variation [T4] . 



In summary, we have produced laser-cooled crystals of 
triply charged 229 Th in a linear Paul trap. Our work 
opens an avenue toward location and precise measure- 
ment of the nuclear transition energy in a single trapped, 
cold ion. Laser excitation of this nuclear transition may 
lead to new levels of precision in timekeeping and tests 
of temporal variation of fundamental constants. 
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Note added. — While writing this manuscript, we be- 
came aware of Ref. [22], where the value of Q = 
3.149±0.032 eb for 229 Th was determined via Coulomb 
excitation of the nucleus, in good agreement with our 
result. 
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